15 openstreetmap

Downloading OpenStreetMap data with a single line of code

Uncomment the following line to install leafmap if needed.

In [1]:

# !pip install geopandas osmnx leafmap

# !pip install geopandas osmnx leafmap

In [2]:

import leafmap

import leafmap

Add OSM data of place(s) by name or ID to the map. Note that the leafmap custom layer control does not support GeoJSON, we need to use the ipyleaflet built-in layer control.

In [3]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_geocode("New York City", layer_name='NYC')
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_geocode("New York City", layer_name='NYC') m

Out[3]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [4]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_geocode("Chicago, Illinois", layer_name='Chicago, IL')
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_geocode("Chicago, Illinois", layer_name='Chicago, IL') m

Out[4]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within boundaries of geocodable place(s) to the map.

In [5]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
place = "Bunker Hill, Los Angeles, California"
tags = {"building": True}
m.add_osm_from_place(place, tags, layer_name="Los Angeles, CA")
m

m = leafmap.Map(toolbar_control=False, layers_control=True) place = "Bunker Hill, Los Angeles, California" tags = {"building": True} m.add_osm_from_place(place, tags, layer_name="Los Angeles, CA") m

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[5]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Show OSM feature tags.

https://wiki.openstreetmap.org/wiki/Map_features

In [6]:

# leafmap.osm_tags_list()

# leafmap.osm_tags_list()

Add OSM entities within some distance N, S, E, W of address to the map.

In [7]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_address(
    address="New York City", 
    tags={"amenity": "bar"}, 
    dist=1500,
    layer_name="NYC bars")
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_address( address="New York City", tags={"amenity": "bar"}, dist=1500, layer_name="NYC bars") m

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[7]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [8]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_address(
    address="New York City", 
    tags={
        "landuse": ["retail", "commercial"],
        "building": True
    },
    dist=1000, layer_name="NYC buildings")
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_address( address="New York City", tags={ "landuse": ["retail", "commercial"], "building": True }, dist=1000, layer_name="NYC buildings") m

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[8]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within a N, S, E, W bounding box to the map.

In [9]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
north, south, east, west = 40.7551, 40.7454, -73.9738, -73.9965
m.add_osm_from_bbox(north, south, east, west, tags={"amenity": "bar"}, layer_name="NYC bars")
m

m = leafmap.Map(toolbar_control=False, layers_control=True) north, south, east, west = 40.7551, 40.7454, -73.9738, -73.9965 m.add_osm_from_bbox(north, south, east, west, tags={"amenity": "bar"}, layer_name="NYC bars") m

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:

2021-11-02 18:51:06 Configured OSMnx 1.1.1
2021-11-02 18:51:06 HTTP response caching is on
2021-11-02 18:51:06 Retrieved response from cache file "cache/8d8366ad6b96fa854a9d6d038c022b8bdcdb6fe6.json"
2021-11-02 18:51:06 Created GeoDataFrame with 1 rows from 1 queries
2021-11-02 18:51:06 Constructed place geometry polygon(s) to query API
2021-11-02 18:51:06 Projected GeoDataFrame to +proj=utm +zone=11 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:51:06 Projected GeoDataFrame to epsg:4326
2021-11-02 18:51:06 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:51:06 Retrieved response from cache file "cache/95fea0c13cd89a78cbaa52ae4b13a8c9ebd9eb21.json"
2021-11-02 18:51:06 Got all geometries data within polygon from API in 1 request(s)
2021-11-02 18:51:06 Converting 1436 elements in JSON responses to geometries
2021-11-02 18:51:06 77 geometries created in the dict
2021-11-02 18:51:06 13 untagged geometries removed
2021-11-02 18:51:06 Created r-tree spatial index for 64 geometries
2021-11-02 18:51:06 Identified 64 geometries inside polygon
2021-11-02 18:51:06 0 geometries removed by the polygon filter
2021-11-02 18:51:06 8 geometries removed by the tag filter
2021-11-02 18:51:06 56 geometries in the final GeoDataFrame
2021-11-02 18:51:07 Pausing 1 seconds before making HTTP GET request
2021-11-02 18:51:08 Get https://nominatim.openstreetmap.org/search?format=json&limit=1&dedupe=0&q=New+York+City with timeout=180
2021-11-02 18:51:08 Downloaded 0.5kB from nominatim.openstreetmap.org
2021-11-02 18:51:08 Saved response to cache file "cache/78258c318356c20c380e1fed844176f90c053cf7.json"
2021-11-02 18:51:08 Geocoded "New York City" to (40.7127281, -74.0060152)
2021-11-02 18:51:08 Created bbox 1500 m from (40.7127281, -74.0060152): 40.72621790503239,40.69923829496761,-73.98821837587228,-74.02381202412771
2021-11-02 18:51:08 Projected GeoDataFrame to +proj=utm +zone=18 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:51:08 Projected GeoDataFrame to epsg:4326
2021-11-02 18:51:08 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:51:08 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:08 Pausing 0 seconds before making HTTP POST request
2021-11-02 18:51:08 Post https://overpass-api.de/api/interpreter?data=%5Bout%3Ajson%5D%5Btimeout%3A180%5D%3B%28%28node%5B%27amenity%27%3D%27bar%27%5D%28poly%3A%2740.699238+-74.023812+40.699238+-73.988218+40.726218+-73.988218+40.726218+-74.023812+40.699238+-74.023812%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27amenity%27%3D%27bar%27%5D%28poly%3A%2740.699238+-74.023812+40.699238+-73.988218+40.726218+-73.988218+40.726218+-74.023812+40.699238+-74.023812%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27amenity%27%3D%27bar%27%5D%28poly%3A%2740.699238+-74.023812+40.699238+-73.988218+40.726218+-73.988218+40.726218+-74.023812+40.699238+-74.023812%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%29%3Bout%3B with timeout=180
2021-11-02 18:51:08 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:10 Downloaded 36.9kB from overpass-api.de
2021-11-02 18:51:10 Saved response to cache file "cache/7d10c2e8d50f39407dda671d56290ea2808b3d4a.json"
2021-11-02 18:51:10 Got all geometries data within polygon from API in 1 request(s)
2021-11-02 18:51:10 Converting 107 elements in JSON responses to geometries
2021-11-02 18:51:10 98 geometries created in the dict
2021-11-02 18:51:10 0 untagged geometries removed
2021-11-02 18:51:10 Created r-tree spatial index for 98 geometries
2021-11-02 18:51:10 Identified 98 geometries inside polygon
2021-11-02 18:51:10 0 geometries removed by the polygon filter
2021-11-02 18:51:10 0 geometries removed by the tag filter
2021-11-02 18:51:10 98 geometries in the final GeoDataFrame
2021-11-02 18:51:10 Retrieved response from cache file "cache/78258c318356c20c380e1fed844176f90c053cf7.json"
2021-11-02 18:51:10 Geocoded "New York City" to (40.7127281, -74.0060152)
2021-11-02 18:51:10 Created bbox 1000 m from (40.7127281, -74.0060152): 40.721721303354926,40.703734896645074,-73.99415065058152,-74.01787974941847
2021-11-02 18:51:10 Projected GeoDataFrame to +proj=utm +zone=18 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:51:10 Projected GeoDataFrame to epsg:4326
2021-11-02 18:51:10 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:51:10 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:11 Pausing 0 seconds before making HTTP POST request
2021-11-02 18:51:11 Post https://overpass-api.de/api/interpreter?data=%5Bout%3Ajson%5D%5Btimeout%3A180%5D%3B%28%28node%5B%27landuse%27%3D%27retail%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27landuse%27%3D%27retail%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27landuse%27%3D%27retail%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28node%5B%27landuse%27%3D%27commercial%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27landuse%27%3D%27commercial%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27landuse%27%3D%27commercial%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28node%5B%27building%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27building%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27building%27%5D%28poly%3A%2740.703735+-74.017880+40.703735+-73.994151+40.721721+-73.994151+40.721721+-74.017880+40.703735+-74.017880%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%29%3Bout%3B with timeout=180
2021-11-02 18:51:11 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:17 Downloaded 2,708.5kB from overpass-api.de
2021-11-02 18:51:17 Saved response to cache file "cache/2c876f231fa35ebe028236a3fb06b51de7794ae0.json"
2021-11-02 18:51:17 Got all geometries data within polygon from API in 1 request(s)
2021-11-02 18:51:17 Converting 20831 elements in JSON responses to geometries
2021-11-02 18:51:17 No outer polygons were created for https://www.openstreetmap.org/relation/3695753
2021-11-02 18:51:17 2922 geometries created in the dict
2021-11-02 18:51:17 70 untagged geometries removed
2021-11-02 18:51:18 Created r-tree spatial index for 2852 geometries
2021-11-02 18:51:18 Identified 2845 geometries inside polygon
2021-11-02 18:51:18 7 geometries removed by the polygon filter
2021-11-02 18:51:18 124 geometries removed by the tag filter
2021-11-02 18:51:18 2727 geometries in the final GeoDataFrame
2021-11-02 18:51:23 Projected GeoDataFrame to +proj=utm +zone=18 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:51:23 Projected GeoDataFrame to epsg:4326
2021-11-02 18:51:23 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:51:23 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:23 Pausing 1 seconds before making HTTP POST request

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[9]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within some distance N, S, E, W of a point to the map.

In [10]:

m = leafmap.Map(center=[46.7808, -96.0156], zoom=12, toolbar_control=False, layers_control=True)
m.add_osm_from_point(center_point=(46.7808, -96.0156), tags={"natural": "water"}, dist=10000, layer_name="Lakes")
m

m = leafmap.Map(center=[46.7808, -96.0156], zoom=12, toolbar_control=False, layers_control=True) m.add_osm_from_point(center_point=(46.7808, -96.0156), tags={"natural": "water"}, dist=10000, layer_name="Lakes") m

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[10]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [11]:

m = leafmap.Map(center=[39.9170, 116.3908], zoom=15, toolbar_control=False, layers_control=True)
m.add_osm_from_point(center_point=(39.9170, 116.3908), tags={"building": True, "natural": "water"}, dist=1000, layer_name="Beijing")
m

m = leafmap.Map(center=[39.9170, 116.3908], zoom=15, toolbar_control=False, layers_control=True) m.add_osm_from_point(center_point=(39.9170, 116.3908), tags={"building": True, "natural": "water"}, dist=1000, layer_name="Beijing") m

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:816: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for merged_inner_linestring in list(merged_inner_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:805: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
  for merged_outer_linestring in list(merged_outer_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/geometries.py:805: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for merged_outer_linestring in list(merged_outer_linestrings):
/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[11]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within the current map view to the map.

In [12]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.set_center(-73.9854, 40.7500, 16)
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.set_center(-73.9854, 40.7500, 16) m

Out[12]:

Make this Notebook Trusted to load map: File -> Trust Notebook

The add_osm_from_view() is only available for the ipyleaflet plotting backend.

In [13]:

# m.add_osm_from_view(tags={"amenity": "bar", "building": True}, layer_name="New York")

# m.add_osm_from_view(tags={"amenity": "bar", "building": True}, layer_name="New York")

Create a GeoPandas GeoDataFrame from place.

In [14]:

gdf = leafmap.osm_gdf_from_place("New York City", tags={"amenity": "bar"})
gdf

gdf = leafmap.osm_gdf_from_place("New York City", tags={"amenity": "bar"}) gdf

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:335: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for polygon in geometry:

2021-11-02 18:51:24 Post https://overpass-api.de/api/interpreter?data=%5Bout%3Ajson%5D%5Btimeout%3A180%5D%3B%28%28node%5B%27amenity%27%3D%27bar%27%5D%28poly%3A%2740.745400+-73.996500+40.745400+-73.973800+40.755100+-73.973800+40.755100+-73.996500+40.745400+-73.996500%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27amenity%27%3D%27bar%27%5D%28poly%3A%2740.745400+-73.996500+40.745400+-73.973800+40.755100+-73.973800+40.755100+-73.996500+40.745400+-73.996500%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27amenity%27%3D%27bar%27%5D%28poly%3A%2740.745400+-73.996500+40.745400+-73.973800+40.755100+-73.973800+40.755100+-73.996500+40.745400+-73.996500%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%29%3Bout%3B with timeout=180
2021-11-02 18:51:24 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:27 Downloaded 16.1kB from overpass-api.de
2021-11-02 18:51:27 Saved response to cache file "cache/fc1b84bb839adfa1ccfe180041a3a3ac4a5fb27b.json"
2021-11-02 18:51:27 Got all geometries data within polygon from API in 1 request(s)
2021-11-02 18:51:27 Converting 49 elements in JSON responses to geometries
2021-11-02 18:51:27 39 geometries created in the dict
2021-11-02 18:51:27 0 untagged geometries removed
2021-11-02 18:51:27 Created r-tree spatial index for 39 geometries
2021-11-02 18:51:27 Identified 39 geometries inside polygon
2021-11-02 18:51:27 0 geometries removed by the polygon filter
2021-11-02 18:51:27 3 geometries removed by the tag filter
2021-11-02 18:51:27 36 geometries in the final GeoDataFrame
2021-11-02 18:51:27 Created bbox 10000 m from (46.7808, -96.0156): 46.87073203354929,46.69086796645071,-95.88427235704086,-96.14692764295916
2021-11-02 18:51:27 Projected GeoDataFrame to +proj=utm +zone=14 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:51:27 Projected GeoDataFrame to epsg:4326
2021-11-02 18:51:27 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:51:27 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:27 Pausing 25 seconds before making HTTP POST request
2021-11-02 18:51:52 Post https://overpass-api.de/api/interpreter?data=%5Bout%3Ajson%5D%5Btimeout%3A180%5D%3B%28%28node%5B%27natural%27%3D%27water%27%5D%28poly%3A%2746.690868+-96.146928+46.690868+-95.884272+46.870732+-95.884272+46.870732+-96.146928+46.690868+-96.146928%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27natural%27%3D%27water%27%5D%28poly%3A%2746.690868+-96.146928+46.690868+-95.884272+46.870732+-95.884272+46.870732+-96.146928+46.690868+-96.146928%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27natural%27%3D%27water%27%5D%28poly%3A%2746.690868+-96.146928+46.690868+-95.884272+46.870732+-95.884272+46.870732+-96.146928+46.690868+-96.146928%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%29%3Bout%3B with timeout=180
2021-11-02 18:51:52 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:54 Downloaded 3,520.6kB from overpass-api.de
2021-11-02 18:51:54 Saved response to cache file "cache/7e1b9b17a87fd207fe2ef4bba35526a26b941cf3.json"
2021-11-02 18:51:54 Got all geometries data within polygon from API in 1 request(s)
2021-11-02 18:51:54 Converting 34091 elements in JSON responses to geometries
2021-11-02 18:51:55 648 geometries created in the dict
2021-11-02 18:51:55 67 untagged geometries removed
2021-11-02 18:51:55 Created r-tree spatial index for 581 geometries
2021-11-02 18:51:55 Identified 581 geometries inside polygon
2021-11-02 18:51:55 0 geometries removed by the polygon filter
2021-11-02 18:51:55 9 geometries removed by the tag filter
2021-11-02 18:51:55 572 geometries in the final GeoDataFrame
2021-11-02 18:51:56 Created bbox 1000 m from (39.917, 116.3908): 39.92599320335493,39.908006796645076,116.40252555265246,116.37907444734753
2021-11-02 18:51:56 Projected GeoDataFrame to +proj=utm +zone=50 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:51:56 Projected GeoDataFrame to epsg:4326
2021-11-02 18:51:56 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:51:56 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:51:56 Pausing 5 seconds before making HTTP POST request
2021-11-02 18:52:01 Post https://overpass-api.de/api/interpreter?data=%5Bout%3Ajson%5D%5Btimeout%3A180%5D%3B%28%28node%5B%27building%27%5D%28poly%3A%2739.908007+116.379074+39.908007+116.402526+39.925993+116.402526+39.925993+116.379074+39.908007+116.379074%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27building%27%5D%28poly%3A%2739.908007+116.379074+39.908007+116.402526+39.925993+116.402526+39.925993+116.379074+39.908007+116.379074%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27building%27%5D%28poly%3A%2739.908007+116.379074+39.908007+116.402526+39.925993+116.402526+39.925993+116.379074+39.908007+116.379074%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28node%5B%27natural%27%3D%27water%27%5D%28poly%3A%2739.908007+116.379074+39.908007+116.402526+39.925993+116.402526+39.925993+116.379074+39.908007+116.379074%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28way%5B%27natural%27%3D%27water%27%5D%28poly%3A%2739.908007+116.379074+39.908007+116.402526+39.925993+116.402526+39.925993+116.379074+39.908007+116.379074%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%28relation%5B%27natural%27%3D%27water%27%5D%28poly%3A%2739.908007+116.379074+39.908007+116.402526+39.925993+116.402526+39.925993+116.379074+39.908007+116.379074%27%29%3B%28._%3B%3E%3B%29%3B%29%3B%29%3Bout%3B with timeout=180
2021-11-02 18:52:01 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:52:04 Downloaded 1,933.1kB from overpass-api.de
2021-11-02 18:52:04 Saved response to cache file "cache/18c0b298eadc540f67bd395df11ce3a66c601ba8.json"
2021-11-02 18:52:04 Got all geometries data within polygon from API in 1 request(s)
2021-11-02 18:52:04 Converting 17009 elements in JSON responses to geometries
2021-11-02 18:52:05 2724 geometries created in the dict
2021-11-02 18:52:05 78 untagged geometries removed
2021-11-02 18:52:05 Created r-tree spatial index for 2646 geometries
2021-11-02 18:52:05 Identified 2645 geometries inside polygon
2021-11-02 18:52:05 1 geometries removed by the polygon filter
2021-11-02 18:52:05 39 geometries removed by the tag filter
2021-11-02 18:52:05 2607 geometries in the final GeoDataFrame
2021-11-02 18:52:07 Configured OSMnx 1.1.1
2021-11-02 18:52:07 HTTP response caching is on
2021-11-02 18:52:07 Retrieved response from cache file "cache/461ac14500db314f1655f9539ed600906dcbe1da.json"
2021-11-02 18:52:07 Created GeoDataFrame with 1 rows from 1 queries
2021-11-02 18:52:07 Constructed place geometry polygon(s) to query API
2021-11-02 18:52:07 Projected GeoDataFrame to +proj=utm +zone=18 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +type=crs
2021-11-02 18:52:07 Projected GeoDataFrame to epsg:4326
2021-11-02 18:52:07 Requesting data within polygon from API in 1 request(s)
2021-11-02 18:52:07 Resolved overpass-api.de to 178.63.48.217
2021-11-02 18:52:08 Pausing 1 seconds before making HTTP POST request

/opt/hostedtoolcache/Python/3.8.12/x64/lib/python3.8/site-packages/osmnx/utils_geo.py:426: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for poly in multipoly:

Out[14]:

		addr:city	addr:country	addr:housenumber	addr:postcode	addr:state	addr:street	amenity	ele	gnis:county_name	gnis:feature_id	...	disused:live_music	disused:amenity	screen	drink:cocktails	building:colour	roof:material	roof:shape	billiards	ways	type
element_type	osmid
node	368043598	New York City	US	53	10014	NY	Christopher Street	bar	7	New York	2082946	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	419366609	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	419367969	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	484078905	Brooklyn	NaN	388	11211	NY	Union Avenue	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	617740843	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
way	664632831	NaN	NaN	1840	10128	NaN	2nd Avenue	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	664910637	Brooklyn	NaN	3	11222	NY	Milton Street	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	665897391	New York	NaN	54	10001	NaN	West 31st Street	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	763525023	NaN	NaN	183	NaN	NaN	East 86th Street	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	16	NaN	NaN
relation	12554650	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	[278963046, 278963089]	multipolygon

1011 rows × 169 columns

---

Last update: 2021-11-02